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Abstract 

Superpartner masses cannot be arbitrarily heavy if super symmetric exten- 
sions of the standard model explain the stability of the gauge hierarchy. 
This ancient and hallowed motivation for weak scale supersymmetry is of- 
ten quoted, yet no reliable determination of this upper limit on superpartner 
masses exists. In this paper we compute upper bounds on superpartner masses 
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in the minimal supersymmetric model, and we identify which values of the 
superpartner masses correspond to the most natural explanation of the hier- 
archy stability. We compare the most natural value of these masses and their 
upper limits to the physics reach of current and future colliders. As a re- 
sult, we find that supersymmetry could explain weak scale stability naturally 
even if no superpartners are discovered at LEP II or the Tevatron (even with 
the Main Injector upgrade). However, we find that supersymmetry cannot 
provide a complete explanation of weak scale stability, if squarks and gluinos 
have masses beyond the physics reach of the LHC. Moreover, in the most nat- 
ural scenarios, many sparticles, for example, charginos, squarks, and gluinos, 
lie within the physics reach of either LEP II or the Tevatron. Our analysis 
determines the most natural value of the chargino (squark) ((gluino)) mass 
consistent with current experimental constraints is ^ 50 (250) ((250)) GeV 
and the corresponding theoretical upper bound is ~ 250 (700) ((800)) GeV. 
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I. INTRODUCTION 



As a candidate for physics beyond the standard model, weak scale supersymmetry has 
several appealing features: It provides an understanding of why a light weak scale is stable, 
it successfully predicts the value of sin^ 6w assuming gauge unification, it predicts a top 
quark Yukawa coupling of order one, leading to a heavy Mt (assuming r lepton and bottom 
quark Yukawa coupling unification), and it provides a natural cold dark matter candidate 
in the form of the lightest superpartner. 

Despite these circumstantial arguments for weak scale supersymmetry, there is not a 
shred of direct experimental evidence to support it. Should we be surprised or discouraged 
that we have not yet found any supersymmetric partners to the standard model particles? To 
date, of the particles we believe to be fundamental, all those observed would be massless if the 
gauge symmetries of the standard model were unbroken. Because the current, experimental 
probes only reach up through the lower fringes of the weak scale, it is not surprising that 
the fundamental particles discovered so far obtain masses as a consequence of spontaneously 
broken gauge symmetries. Their superpartners, by contrast, can have gauge invariant mass 
terms, provided supersymmetry is broken. Although they are not required to be light by 
gauge symmetries, there is a theoretical upper limit on their masses above which the weak 
scale does not arise naturally. As the scale of supersymmetry breaking is increased, the weak 
scale can only remain light by virtue of an increasingly delicate cancellation. Requiring that 
the weak scale arises naturally places an upper bound on superpartner masses. 

In this paper, we attempt to quantify the relationship between naturalness and superpart- 
ner masses. Using recently formulated naturalness measures we compute the most natural 
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value of the superpartner masses, the extent to which naturalness is lost as experimental 
bounds on superpartner masses increase, and a theoretical upper limit to the masses of 
superparticles. 

In section two we review the naturalness measures used in our study. Section three 
is devoted to a discussion of radiative electroweak symmetry breaking in the minimal su- 
persymmetric extension of the Standard Model (MSSM) and to details of the numerical 
methods we employed in our analysis. The results, presented in section four, demonstrate 
that the MSSM can not accommodate the weak scale naturally if superpartner masses he 
beyond the reach of LHC. Moreover, in the most natural cases, physics beyond the standard 
model has a good chance of being discovered at LEP II or the Tevatron. 



II. MEASURING FINE TUNING 

In this section we review the recently formulated naturalness measures we use in our anal- 
ysis. A more detailed motivation and derivation of these criteria can be found in Ref. [1]. 
Any measure of naturalness contains assumptions about how the fundamental parameters of 
a Lagrangian are distributed. If we parameterize these assumptions, a quantitative measure 
of naturalness follows directly. Consider a Lagrangian density written in terms of fundamen- 
tal couplings specified at the high energy boundary of the effective theory: £(01,02, ■■■an)- 
At a low energy scale, we can write the Lagrangian in terms of physical observables X {e.g., 
X — M|). These observables will depend on the Oj through the renormalization group 
and possibly on a set of minimization conditions: X — X{a)^ If we assume the probability 
distribution of a fundamental Lagrangian parameter a is given by 

,„/ X f(a)da 

dP{a) = )\ , 2.1 
S f{a)da 

a likelihood distribution for the low energy observable X follows 

/ f{a)da = / p{X)dX . (2.2) 

J a- J X{a-) 
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The value of an observable X is unnatural if it is relatively unlikely to end up in an interval 
u{X) about X compared to similarly defined intervals around other values of X. The 
probability that X lies within an interval u{X) about X has weight up. So we define our 
quantitative measure of naturalness as 

{up) 

^ = ^(xM^y (2.3) 

where 

The conventional sense of naturalness for hierarchy problems corresponds to an interval 
■u = X. 1^ With this prescription, fine tuning corresponds to 7 >> 1. The 7 defined in Eq. 
(2.3) is proportional to the Barbieri-Giudice sensitivity parameter c{X, a) = \ {a/X){dX/da)\ 
0). We can use Eqs. (2.3-4) to define an average sensitivity c through the relation 

7 = c/c . (2.5) 

This definition of c gives 

J daaf{a)c{X;a)-^ 

i/c = 7rm • ^^-^^ 

a J [a) J da 

The naturalness measures defined by Eqs. (2.5-6) are a refinement of Susskind's description 
of Wilson's naturalness criteria 0: Observable properties of a system, i.e., X, should not 
be unusually unstable with respect to minute variations in the fundamental parameters, a. 
In other words, X{a) is fine tuned if the values of the fundamental parameters a are chosen 



§ For example, in a theory of fundamental scalars, the scalar mass is related to the cut-off A and 
the bare term niQ by: = g'^A'^ ~^o- this theory we must adjust with the same precision to 
place the scalar mass squared in a 1 (GeV)^ window whether the scalar mass is 0(A) or O(10~^^A). 
A small mass for the scalar is unnatural in the sense that a small change in leads to a large 
fractional change in . 



so that X depends on the a in an unusually sensitive manner when compared to other values 
of the fundamental parameters a. Sensitivity in this case is understood to mean that a small 
fractional change in a leads to a large fractional change in X. Q 

Returning to Eq. (2.4-2.6), we see that three choices need to be specified before we 
can make practical use of this prescription. First, the choice of /(a) reflects our theoretical 
prejudice about what constitutes a natural value of the Lagrangian parameter a. We will 
make two different choices for /(a) as an aid in determining how sensitively the bounds 
we derive depend on this theoretical prejudice: /(a) = 1 and /(a) = 1/a. We denote 
the corresponding naturalness measures by 71 and 72, respectively. The bounds we derive 
on superpartner masses in section four are fairly insensitive to this choice. Second, the 
conventional notion of naturalness for hierarchy problems is u{X) = X This choice has 
already been made in Eq. (2.6) and it is implicit in the qualitative statement of naturalness 
written above. Finally, the range of integration (a_,a+) for the averaging must be chosen. 
This range will be discussed in section four. 

III. THE MSSM 

All the chiral interactions of the MSSM are described by its superpotential 

W = TIYuKQ + ^Yd^dQ + tYe<lrfL + fiK^d . (3.1) 

The /x-term explicitly breaks the Peccei-Quinn symmetry and avoids a phenomenologically 
disastrous axion. In addition to all the particles of the SM, there are thirty-one new ones 
including three new Higgs bosons. 

Supersymmetry is explicitly broken in the MSSM using soft terms derived from the low 
energy limit of supergravity (SUGRA) theory. The form of the soft SUSY breaking potential 

** In deriving the naturalness criteria Eqs (2.3-2.6), we have attempted to make explicit the discre- 
tionary choices inherent in any quantitative measure of naturalness. In any particular application, 
in order to obtain a reliable measure of naturalness, these choices must be made sensibly. 
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in the MSSM includes mass terms for all the scalars and for the gauginos as well as bilinear 
and trilinear terms following from the Kahler potential of the SUGRA theory in the low 
energy limit. 

A generic feature of minimal low energy SUGRA models is universality of the soft terms. 
Universality implies that all the scalar mass parameters are equal to the gravitino mass, 
mo, at some high energy scale which we take to be the scale of gauge coupling unification, 
Mx = 10^^ GeV. All soft trilinear couplings share a common value, Aq, that can be related 
to the soft bilinear coupling, Sq, depending on the form of the Kahler potential. To some 
extent, universality in the soft breaking terms is required in order to avoid unwanted flavor 
changing neutral current effects. Since the gauge couplings unify, the gaugino mass parame- 
ters are assumed equal to a common value, mi/2, at Mx- Consequently, the minimal model 
introduces five new parameters, mg, Aq, "^1/2, -Bo; and /iq. However, it is very predictive 
since these account for the masses of thirty-one new particles [Q. 

In the MSSM, the electroweak symmetry is broken radiatively [§-§|. In our analysis, we 
use the 1-loop effective Higgs potential 

^i-ioop = K) + AV^i , (3.2) 
where the expression for the 1-loop correction is given by 

A^i = ^E(-l)"'(2^^ + lK(lngf-2) ■ (2.3) 

The rrii represent the field dependent masses of the particles of the model and the Sj the 
associated spins. We include the contributions of all the MSSM particles in the 1-loop 
correction. 

Using the renormalization group, the parameters are evolved to low energies where the 
potential attains validity. This RG improvement uncovers electroweak symmetry breaking. 
The exact low energy scale at which to minimize is unimportant as long as the 1-loop effective 
potential is used and the scale is in the expected electroweak range. The minimization scale 
will arbitrarily be taken to be Mz- If the electroweak symmetry is broken, minimization 
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yields non-zero values for the vacuum expectation values (VEVs) of the two Higgs fields, v^i 



and Vd, or equivalently v = ^v'^ + and tan/3 = Vu/vd- The two minimization conditions 
may be expressed as follows 



^\M,) = "''\S. : - >^ ' (3-4) 



mi — mi tan^ 3 1 

tan2/5-l 2' 

(ml + mL + 2/i2) sin 2(3 , ^ 

B(Mz) = ^7 ^ , 3.5 

where ^ = ^ + dAVi/dv^ ^. In the process of integrating the 2-loop renormalization 
group equations, the threshold corrections due to all the light particles are implemented as 
step functions 0. 

The procedure we follow to analyze the MSSM assumes the following 4+1 free input 
parameters: Aq, mo, mi/2, tan/3o, and sign(/i) since it is undetermined from Eq. (3.4). 
The other parameters of the MSSM are fixed using the following constraints: Electroweak 
breaking in the form of two minimization conditions at M^, the physical masses of the 
bottom quark and r lepton, and the value of the strong coupling at Mz- Therefore, solutions 
for Bq, fiQ, yT-{Mx), a3{Mx), and Mf are found consistent with the RG, the above constraints, 
and specified values for the free input parameters. We take the value of the strong coupling 
at Mz to be a^i^Mz) = .118. The corresponding value of the strong coupling at Mx is 
determined based on this low energy constraint. The values of ai{Mx) and a2{Mx) are 
set equal and fixed at 1/25.3. This constant value for ai^2 at Mx never leads to more than 
about 1% and 3% error in ctem and siii^Ow, respectively. The difference in a^i^Mx) and 
C(i,2{Mx) is at most 3% and can be accommodated using GUT thresholds. 

Not all input values for the free parameters will yield adequate solutions, and the 4 + 1 
dimensional parameter space must be explored and restricted using various criteria. Cases 
are rejected based on the existence of color/charge breaking vacua or a charged lightest 
supersymmetric particle (LSP). In arriving at the superpartner mass bounds, the fine tuning 
prescription, Eq. (2.3), is applied to all solutions found in a grid of approximately 2000 
points bounded as follows: \Ao\ < 400 GeV, < mo < 400 GeV, |mi/2| < 500 GeV, 
1 < tan/5(Mx) < 15, and sign(/i) = ±. 
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IV. ANALYSIS 



The essential, novel feature of the fine tuning measure 7 is to evaluate the sensitivity, 
c, of a physical quantity relative to a benchmark, c. We have derived a formula for this 
benchmark in section three and in Ref. [1]. This prescription for calculating c requires us 
to choose a range of integration (a_, a+). We use two conditions to define a suitable range 
of integration. First, we integrate over the all values of a where SU{'i) x SU{2) x f/(l) is 
broken to SU{'i) x U{l)em- The resulting limits on the range of integration generally come 
from two conditions on the value of Mz- The minimum value of Mz cannot be less than 0, 
and its maximum value cannot exceed some upper bound, often set by the requirement that 
sneutrino squared masses be positive. Second, in our analysis we only consider points where 
we are able to find a significantly large range of integration. If the range of integration is 
not suitably large, we will fail in our attempt to compare the sensitivity of Mz-, when a is 
chosen so that the value of Mz is 91.2 GeV, to the average sensitivity. Inspection of Eq. 
(2.6) shows that in the limit of vanishing (a+ — a_), c approaches c, and 7 tends to one. To 
eliminate spurious calculations of 7, we only consider cases were 5a = a+ — a_ exceeds a/4 
or a/8 for M^(a) = 91.2 GeV. We find that typically this has the effect of removing points 
where S'[/(3) only remains unbroken as the result of a fine tuning. 

Figures 2-9 display correlations between the superpartner masses and fine tuning. For 
each solution point, we computed the fine tuning with respect to the common scalar mass, 
the top quark Yukawa coupling, and the common gaugino mass. Then, for each individual 
solution, we define 7 as the largest of these fine tunings. Many earlier studies of naturalness, 
as well as employing measures of sensitivity instead of fine-tuning, considered the naturalness 
of the Z mass with respect to individual parameters separately. This separation can lead 
to a significant underestimate of fine tuning. In particular, we have compared the lower 

''^"'"In fact, the original bound of c < 10 imposed by Barbieri and Giudice can no longer be satisfied. 
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envelopes defined by scatter plots, and we find explicitly that, if fine tuning is plotted as 
a function of a particular coupling or mass, the envelope defined by 7 cannot in general 
be constructed from the individual envelopes for 7(mo), 1(^1/2), and 7(?/t). Figures 
2-9 display the fine tuning measure 7 plotted against selected superpartner masses. The 
individual points shown in these figures correspond to the grid of approximately 2000 points 
discussed in section three. We caution the reader that the density of these points is not 
an indication of how likely particular values of the superpartner masses or 7 are. This is 
because the grid we have used is not completely uniform, and more importantly because 
the minimization conditions (3.4-5) have been used to determine the values of Bq and fiQ. 
The dashed and dotted curves in these figures show the minimum fine tuning necessary for 
a particular value of the superpartner mass. The likelihood or naturalness of a particular 
value of a superpartner mass scales like I/7. 

Figure 1 contrasts the sensitivity parameter c with our measure of fine tuning 7. We 
see that currently viable solutions depend on at least one fundamental parameter in a fairly 
sensitive manner, however the fine tuning curve, 7, shows that this sensitivity is not always 
unusual. 

Figures 2a-2b display the correlation between the gluino mass and the fine tuning param- 
eters 7i and 72. This plot and, unless otherwise noted, the following plots are constructed 
from solution points consistent with the current LEP limits on superpartner masses |]lOl . We 
have taken the limits on the sneutrino and the charged superpartner masses to be M^/2, 
and the lower limit on the light Higgs mass as 60 GeV. If no superpartner masses lie below 



A calculation of the sensitivity of tlie Z mass with respect to mo, mi/2, Ut, and 53 gives c > 30 
(see Fig. 1). 

is a reflection of the fact that because the Z mass depends on several parameters, even if 
another variable is fixed, it is easy to find solutions where the Z boson's dependence on an isolated 
fundamental parameter is relatively insensitive. 



8 



these limits the most natural value of the gluino mass is about 260 GeV, above the pub- 
lished CDF limit of 141 GeV and also above the recently reported limit from D0 |Tl|. For 
potential, future search limits at the Tevatron see for example Ref. [12]. If we require that 
fine tunings are at most a part in ten, the gluino mass should lie below ~ 600 — 800 GeV, a 
value that should be easily accessible at the LHC [IT3 . 



Figures 3a-b displays the correlation between fine tuning and the lightest squark mass of 
the first and second generation. The analogous plots for the top squark mass are shown in 
Figs. 4a-b. The most natural value of the stop mass is around 220 GeV and for the lightest 
of the remaining squarks it is about 240 GeV. This is close to the preliminary mass limit 



reported by D0 at Glasgow ||rT|. If we require that fine tunings are at most a part in ten, 
the stop mass should lie below ~ 500 — 600 GeV and the lightest of the remaining squark 
masses should lie below ~ 600 — 800 GeV. 

Figures 5a-b display the correlation between fine tuning and the lightest chargino mass. 
This plot displays solution points consistent with the LEP derived constraints on superpart- 
ner masses with the exception of the chargino mass. The most natural value of the lightest 
chargino mass, corresponding to the smallest 7, is around 50 GeV. Note that a significant 
region of the most natural solutions lie within the physics reach of LEP II, which should be 
able to search for charged particles up to the kinematic limit [|14|. The lightest chargino 



mass should not exceed ~ 200 — 300 GeV if 7 < 10. 

Figures 6a-b display the correlation between fine tuning and the mass of the lightest 
superpartner. The most natural value of the LSP mass appears to be around 42 GeV, and 
the theoretically favored values of the LSP mass are concentrated below 70 GeV. The LSP 
can not be heavier than 150 GeV if 7 < 10. This bound provides a more stringent limit 
than bounds set by the requirement that the LSP not over-close the universe. 

Figures 7a-b summarizes the mass predictions for all the superpartners. The upper and 
lower ends of the bars correspond to 7 < 10 and the current experimental limits, respectively. 
The diamond point represents the 7 < 5 mass limit, and the square represents the most 
natural value for the respective sparticle mass. 



Finally, for completeness we display the correlation between the lower bound on fine 
tuning and the fundamental parameters mo and |/io| in Figs. 8 and 9. 



V. CONCLUSIONS 

As the mass limits on superpartners increase, it becomes increasingly difficult to accom- 
modate a light weak scale naturally. We have presented a detailed study of the relationship 
between superpartner masses and naturalness. This analysis demonstrates that supersym- 
metry can not accommodate the weak scale without significant fine tuning if superpartner 
masses lie beyond the physics reach of the LHC In addition our analysis reveals that the 
most natural values of these masses often lie well below 1 TeV. We note that our limits are 
higher than those which would be obtained using conventional sensitivity criteria, but they 
lie below the bounds found in common folklore. In light of our results, we feel the poten- 
tial for the discovery of physics beyond the standard model before the LHC is promising. 
However, this optimism should not be interpreted as a guarantee that LEP II or the Teva- 
tron will see superpartners even in the case when supersymmetry is relevant to electroweak 
symmetry breaking. A more detailed application of these naturalness measures to collider 
SUSY discovery reaches is in progress ll^ . 
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Figure 1: Curves representing the lower envelope of regions defined by 

c = max{c(mi/2), c{mQ),c{yt), c{gs)} and 

7i,2 = max{7i,2(mi/2),7i,2c(mo),7i,2(?/t),7i,2(5'3)} plotted as a function of tan/3. The 
upper curve represents the amount of sensitivity required by current experimental 
superpartner limits, and the lower curves display the amount of fine tuning. 

Figures 2a-b: The fine tuning measures 71^2 as a function of the gluino mass. 

Figures 3a-b: The fine tuning measures 71^2 as a function of the lightest squark mass of the 
first two generations. 

Figures 40'-b: The fine tuning measures 71^2 as a function of the lightest top squark mass. 

Figures 5a-b: The fine tuning measures 71^2 as a function of the lightest chargino mass. 

Figures 6a-b: The fine tuning measures 71^2 as a function of the lightest sparticle mass. 

Figures 7a-b: Superpartner mass ranges. The upper and lower ends of the bars correspond 
to 7 < 10 and the current experimental limits, respectively. The diamond (square) 
represents the limit 7 < 5 (the most natural value). 

Figures 8a-b: The fine tuning measures 71,2 as a function of the common scalar mass mo. 

Figures 9a-b: The fine tuning measures 71^2 as a function of the mixing parameter |/io|- 
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